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O \ Abstract 

(N 

We investigate the totally asymmetric simple exclusion process (TASEP) in the 

\^ \ presence of a bottleneck, i.e. a sequence of consecutive defect sites with reduced 

W ■ hopping rate. The influence of such a bottleneck on the phase diagram is studied by 

^ . computer simulations and a novel analytical approach. We flnd a clear dependence 

of the current and the properties of the phase diagram not only on the length of 

the bottleneck, but also on its position. For bottlenecks near the boundaries, this 

motivates the concept of effective boundary rates. Furthermore the inclusion of a 

second, smaller bottleneck far from the flrst one has no influence on the transport 

(^T) I capacity. These results will form the basis of an effective description of the disor- 

O ■ dered TASEP and are relevant for the modelling of protein synthesis or intracellular 

transport systems where the motion of molecular motors is hindered by immobile 

blocking molecules. 
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1 Introduction 



Driven diffusive systems play an important role in statistical physics. They 
not only serve as paradigm for non-equilibrium behaviour [UEEJ, but also as 
models for transport processes like vehicular traflSc [3], granular flow through 
narrow pipes [5j and biological transport by motor proteins PITEEj. The 
simplest of these models is the totally asymmetric simple exclusion process 
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(TASEP) which was first introduced to describe protein polymerization in ri- 
bosomes [TO] (for a recent review, see [H]). The TASEP exhibits some generic 
properties hke boundary induced phase transitions [il2j that also occur in more 
complex driven systems. It was solved exactly p^SlfTl] and results can be used 
for qualitative and quantitative approaches to other systems. The phase dia- 
gram consists of three phases (high density (h), low density (1), and maximum 
current phase (m)), depending on whether the current is limited by the par- 
ticle input, output or the transport capacity of the bulk, and is generic for a 
large class of driven diffusive systems [i5\ . 

While homogeneous systems are extensively investigated, systems with inho- 
mogeneous transition rates are not that well understood, especially for the 
case of open boundary conditions. It is well-known that in nonequilibrium 
systems even weak disorder can lead to drastic changes [T6||T7j. For driven 
diffusive systems already a single defect site^ i.e. a site with reduced hopping 
rate, can have a global effect on the stationary state, see e.g. p^81[T9] for the 
case of the TASEP with periodic boundary conditions. Here it was shown that 
in an intermediate density regime \p— 1/2 1 < Pc the current becomes density- 
independent, i.e. the flow-density relation (fundamental diagram) exhibits a 
plateau. In contrast, in the limits of low and high densities, the current is 
unchanged by the presence of the defect. 

Here we consider the TASEP with open boundaries in the presence of a bot- 
tleneck of length /, i.e. a sequence of / consecutive "slow" bondilJ, where 
the hopping rate of the particles is reduced compared to the rest of the sys- 
tem. We will study the effect of such bottlenecks on the phase diagram of 
the TASEP, focussing on the maximal current that can be maintained by the 
system through optimization of the input and output rates, i.e. its transport 
capacity. 

In open systems, the main effect of the bottleneck (or defects in general) is a 
decrease of the transport capacity [20j . This leads to an enlarged maximum 
current phase compared to the pure system. In addition the density proflles in 
the maximum current phase exhibit phase separation into a high density and 
low density regime separated by a shock. 

The effect of defects in systems with open boundaries has been studied earlier 
by several authors. So far no exact solutions are known and one has to rely on 
computer simulations with Monte Carlo (MC) methods and mean-fleld type 
approximations. 

Kolomeisky [20] investigated the TASEP in the presence of a single defect 
site (/ = 1) deep in the bulk. Analytical results are obtained by a mean fleld 
approach which neglects correlations on the slow bond by dividing the system 



A bottleneck of length / = 1 corresponds to a single defect. 



into two homogeneous ones coupled at the defect site. Therefore this approach 
is hmited to single defect sites (/ = 1). Comparison with MC results shows a 
good agreement, at least for the high- and low-density phases. The mean field 
treatment also leads to deviations of density profiles near the defect site from 
the MC results. 

In [21] , Chou and Lakatos gave an analytical approach treating a finite number 
of bottlenecks of arbitrary lengths. Their finite segment mean field theory 
(FSMFT) divides the system into segments which contain one or a few defect 
sites. By determining the leading eigenvector of the corresponding transition 
matrix and matching the currents of the diflFerent subsystems, predictions for 
the current through the system are obtained which at the same time take 
into account correlations near the bottlenecks. Although the method allows 
to treat several defect sites of arbitrary hopping rates within the segment, 
it is restricted to segment lengths of less than 20 sites due to the numerical 
complexity. 

The results for one and two defect sites (/ = 1) in the bulk were recently 
extended by Dong et al. [22]. Additionally they have considered the case where 
a single defect is located near the boundary of the system and found an edge 
effect induced by the interaction of this defect with the boundary. This implies 
a dependence of the current on the position of the defect. Finally, the effects 
of Langmuir kinetics, i.e. particle-creation and -annihilation in the bulk have 
been shown [23] to lead to a rich phase diagram including novel phases. 



The present investigation generalizes the previous works in several aspects. 
Especially we will study systematically the dependence of the current on the 
length / of the bottleneck and its position. For this purpose we develop an 
analytical approach to calculate the current and critical boundary rates, called 
the interacting subsystem approximation (ISA). Mathematically it requires 
finding a specific root, constrained by physical requirements, of a polynomial 
which has a degree of order /. For longer bottlenecks this approach is more 
eflBcient than the approach of [21] which demands to search an eigenvector of a 
2^ X 2^-matrix. This allows us to generalize the single defect results of Dong et 
al. [22] to longer bottlenecks. The increase of the maximum current by moving 
defects near the boundary is reproduced by ISA and extended to the case / > 1 
[positive edge effect). However, for lower entry rates ISA predicts a decrease 
of the current [negative edge effect) which is confirmed by MC simulations. 
We also observe in MC simulations an intermediate regime which shows a 
non-monotonic dependence of the current on the position of the bottleneck. 
Furthermore we conclude that for bottlenecks near the boundaries there is no 
phase transition to a maximum current phase, since the current is not constant 
in any regime, but it approaches the maximum current asymptotically for 
higher entry and exit rates, and there is no macroscopic phase separation. 



We also consider the case of two bottlenecks of arbitrary lengths and sep- 
aration generalizing the corresponding results of [22] for two single defects. 
Therefore we also investigate the case that one of the bottlenecks is near the 
boundaries. The results will motivate a concept of effective boundary rates 
that encompasses the effect of boundary-near bottlenecks. This provides the 
basis for an extensive study [21] of finite defect densities, i.e. a macroscopic 
number of slow bonds. Here previous investigations for periodic [25ll26ll27l[25] 
and open systems [2 8l29ll3QII31ll32|[33 ] have revealed surprising results. For the 
open system, for example, it has been found [30j that the position of the phase 
transitions is sensitively sample-dependent even for large systems. 



2 Definition of the model 



We consider a TASEP consisting of L sites which can either be empty or 
occupied by one particle. Throughout this paper we are mainly interested in 
the steady- state in the limit L ^ oc. With each bond between neighbouring 
sites j and j + 1 we associate a hopping rate Pjj+i which corresponds to the 
rate at which a particle will move from j to its right neighbour j + 1 if this is 
emptjo. At the boundary sites j = I and j = L particles can be inserted and 
removed, respectively. If site 1 is empty a particle will be inserted there with 
rate a. On the other hand, if site L is occupied this particle will be removed 
with rate /?. Here we will use a random-sequential update corresponding to 
continuous-time dynamics. 

In this paper we study the eflFect of a single bottleneck of length /, i.e. a section 
of/ consecutive slow bonds (or sites) with hopping ratepjj+i = q < 1 (Fig.d]). 
All other bonds are fast bonds with hopping rate Pjj+i = 1. 

In order to determine relevant quantities like the current and critical boundary 
rates, we approach the problem by virtually dividing the system into three 
interacting subsystems (Fig. [1]): A system of Li sites with hopping rate pj = 1 
at the left end of the chain, a section of length L2 = / + 1 consisting of the 
bottleneck with / slow sites (hopping rate q) extending from site s = Li + 1 
to site Li + I = Li + L2 — I plus an extra fast site / = Li + L2 at the right 
end, and finally another pure system of length L3 with hopping rates pj = 1. 
Note that subsystem 2, which contains the bottleneck, starts and ends with 
a fast bond. This symmetry of the subsystem will be more convenient for our 
investigations. 



^ SorQetimes it is more convenient to associate hopping rates pj with sites j. These 
site rates are related to the bond rates by pj = Pj,j+i- We will use both terminologies 
here. 
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Fig. 1. Illustration of the TASEP with a bottleneck of length /, i.e. / consecutive 
defect bonds (or sites resp.). The defect bonds/sites are characterized by a hopping 
rate q < I, whereas at the other bonds/sites the hopping rate is p = 1. The system is 
virtually divided into three subsystems with virtual entry/exit rates corresponding 
to the mean occupation/ vacancy on the adjacent sites. 

The main idea of our approach is to use the exact solution for the station- 
ary state of the pure TASEP with random-sequential dynamics p!3l[T4] . This 
solution applies to all three pure subsystems defined above. The interactions 
between these subsystems are described by suitably chosen boundary rates. 
We neglect correlations at the sites connecting the subsystems. Introducing 
the occupation number Tj = 0, 1 of site j, this explicitly means that we assume 



(^j(l - ^j+i)) - Pj(l - Pj+i) fOT J = s-lj 



(1) 



where pj := (tj) is the average local density at site j. This approach is similar 
in spirit to that of Kolomeisky [20] for the case of a single defect (/ = !). 
However, there correlations on the slow bond are neglected whereas here we 
will treat the bottleneck also as a pure system (of reduced hopping rate) so 
that the most relevant correlations induced by it are taken into account. The 
three subsystems are coupled through virtual boundary rates aj , /3j . Note that 
all these rates are associated with fast bonds. In the following we will call this 
approach the interacting subsystem approximation (ISA) to distinguish it from 
the usual mean- field approach which neglects all correlations. 

The rate equations at the sites connecting the subsystems are then given by 



■^{ts) = ^2(1 - (ts)) - q{rs{l - Ts+i)) , 



(2) 
(3) 



and 



|(r;)=g(r;_i(l-r;))-/?2(r;), (4) 

— (r/+i) =0^3(1 - (r/+i)) - (r/+i(l - r/+2)) , (5) 

with the virtual boundary rates 



012 = Ps-i , as = pf , (3i = l- ps, /?2 = 1 - P/+1 , (6) 

in terms of the average local density pj. For completeness, we also define ai := 
a and /Ja := /?. Note that the mean-field factorization of expectation values 
only applies to the sites at the boundaries of the subsystems. All correlations 
within the subsystems will be taken into account exactly! 

From the exact solution of the TASEP with A^ sites for hopping rate p = 1, 
we know the "partition function" ^13j 



which is exact for all system sizes A > 1. Here V{{Ti}) is the probability 
of finding the stationary system in a configuration {r^}. The singularity at 
a = f3 can be removed by taking the limit a ^ f3 which yields the analytic 
continuation 



Z(a, N) := Ito Z(a; 0, N) = g '-%^^U + D Q' ■ (8) 



The corresponding results for the TASEP with hopping rate p ^ I can be 
obtained by rescaling of the boundary rates a -^ a/p^ (3 -^ (3 /p. Thus the 
partition function in a pure TASEP with hopping rate p is 

Z,{a^(5,N) = z(-A.N] . (9) 

\P P J 



This result will be used in the following for different values of A, depending 
on the length of the subsystems. 

The exact current Jq of a pure system of length A, which is the quantity most 
relevant in the following, can be expressed through the partition function 
as [13] 

Jo(a.A/V) = ^'°;f-f-" . (10) 



Conservation of current then yields the central equations of the ISA 

Jo(q^i, A, ^i) = -^0 — , — , ^2 = -^o(<^3, /?3, ^3) (11) 

\q q J 



which express the fact that the current in the stationary state is the same 
in all three interacting subsystems. Note that we have used ([9]) and f ITOl ) to 
express the current in the bottleneck (subsystem 2) through the result Jq for 
a pure system. 

Inserting ([7j) and f ITOl ) in f ITTI ) we obtain a set of two algebraic equations with 
six variables. In principle one can solve these equations, if four of the variables 
are known or more equations are given that determine variables. 



3 One bottleneck far from the boundaries 



First we study the case where the bottleneck is far from the boundaries. So 
we assume that Li and L3 are large. From the latter condition it is expected 
that the topology of the phase diagram is the same as that of a system with 
a single defect (/ = !). We can therefore follow [20j and classify the phases by 
the phases of the pure subsystems 1 and 3 which are 0{L). 

Though at first glance one could expect nine possible phases corresponding 
to all possible combinations of three phases l,h,m (see Sec. [1]) that can be 
realized in the pure subsystems 1 and 3, it was argued in [20j that only the 
combinations 1-1, h-h and h-1 can exist. The 1-1-phase corresponds to a global 
low density phase (L)^ while h-h corresponds to a high density phase (H). In 
both cases, only around the bottleneck there are local deviations from the 
density profile of the homogeneous system. In the h-1-phase, we have phase 
separation which cannot be observed in the pure system. The current in the h- 
phase is independent of the entry rate, while in the 1-phase it is independent of 
the exit rate. Thus in a phase separated h- 1-phase, the current is independent of 
both boundary conditions and takes a maximum value (see below). Although 
it is sometimes called a maximum current phase (M) like in the pure system, 
its properties differ from the maximum current phase in the pure system not 
only by the occurrence of phase separation, but also by the absence of algebraic 
boundary layers. Therefore we prefer the terminology phase separated regime 
(PS). Furthermore, the transition to this phase corresponds to a transition of 
subsystem 1 from 1 to h and vice versa respectively, which is accompanied by 
a discontinuity of the mean density (p) = X^j Pj- According to this, it can be 
classified as a first order transition in contrast to the pure system where the 
transition to the maximum current phase is of second order. 



Though the system with a bottleneck is not exactly invariant under the particle- 
hole symmetry operation defined by Tj ^^ 1 — Tj^ a ^^ (3^ j ^^ L — j^ this trans- 
formation only changes the position of the bottleneck, but still leaves it far 
from the boundary. Thus the particle-hole transformation leaves the phases 
of the subsystems unchanged, since it only changes their sizes, but they stay 
0{L). Therefore, we can conclude that the phase diagram must be symmet- 
ric with respect to the line a = jS. This symmetry constraint yields that the 
transition line between high and low density phase must be at a = /?. 

We now want to determine the critical entry rate a* at the transition from the 
L-phase to the PS-phase for fixed (3 in terms of the analytical ISA approach 
introduced in the last section. It corresponds to the transition of (the pure) 
subsystem 1 from low to high density phase which occurs at a* = ai = /?i. 
Furthermore we know that at this point J = ai{l — ai). Subsystem 3 still 
remains in the low density phase. Therefore we also have J = 0^3(1 — 0^3). Since 
the current must be the same as in subystem 1 and a^ must be smaller than 
1/2 we conclude 0^3 = a*. From the definition of the virtual boundary rates 
(ED, (E]) we obtain by a simple transformation 



f32 = J/pf = J/a'' = l-a\ (12) 

0^2 = J/ {I - Ps) = J/ Pi = J/o^" = l-a\ (13) 

where the first equality in each equation can be found e.g. in [13]. From the 
exact solution of the pure TASEP (ITOl) . ([71) and the conservation of current 
(HH), it follows 



/1-a* 1-a* \ 

J , , / + 1 = a*(l - a*) 

\ q q J 

^^ a*(l - a*)Z((l - a*)/g, / + 1) = Z((l - a*)/g, /) , (14) 

where Z{1 — a, / + 1) is the limit defined in (Ej) and Z{1 — a, / + 1) 7^ for 
< a < 1. 

Eq. ( iT4l) is essentially a polynomial in 1/a* (or a*, respectively) and can be 
solved numerically (or analytically for small values of /). The requirement 
that a physical relevant solution for a* has to be in the interval [0, |] gives a 
unique solution Qf*(g, /) for the transition point to the PS-phase. In general, 
the solution depends on the length / of the bottleneck and the slow hopping 
rate q. For / = 1, for example, equation (IT41) can be transformed into 

2a*^ - (2 + 3g)a* + 2? = . (15) 



The relevant solution is then 
2 + 3g /(2 + 3g)2 



a 



16 



<?• 



(16) 



Explicitly, for the value q = 0.6 used in most simulations, one obtains tt* = | 
by evaluating (IT61) . 

The current at the transition point corresponds to the maximum current J* in 
the PS-phase which can be interpreted as the transport capacity of the system. 
From (IT^ we can see that 



r{q,l) = a*{qj){l-a*{qj)), 



(17) 



which yields J* (0.6, 1) = ^ for q = 0.6 after taking into account (IT6l) . 

In Fig. [2] we have plotted the dependence of J*(g, /) on the slow hopping rate 
q (for fixed /) obtained by (IT7I) and compare these analytical ISA-results with 
results from Monte Carlo (MC) simulations for different bottleneck lengths. 
For a single defect site (/ = 1) we have compared the analytical results obtained 
by (1161) and (IT7I) with the results obtained by pure mean field approximation 
Jmf = (i-^ )2 [20] . Obviously also for / = 1 the results obtained by f lTTI ) 
are more accurate than mean field results, while for longer bottlenecks, to 
our knowledge, no proper mean field approximations are known. Note that 
our approximation takes into account correlations on the slow bonds, only 
correlations on sites adiacent to the bottleneck are neglected. 
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Fig. 2. Maximum current in a system with a bottleneck far from the boundaries in 
dependence on the slow hopping rate q. It is determined by Monte Carlo simulations 
for a = /? = 0.5 and by solving the ISA equation (fMll for a bottleneck of length / = 1 
(left) and / = 6 (right). For comparison the mean field approximation J^p = .-^^ y2 
for / = 1 is included. 



Fig. [3] shows the maximum current in dependence of the bottleneck length /. 
Here we compare MC-simulations with the results obtained by f lTTI ) . The ISA 



results systematically underestimate the real current. The deviation is largest 
for small /, but it does not exceed 3%. For larger bottlenecks, the agreement 
improves. We also note that we have chosen q = 0.6 since for this value the 
observed deviations have been found to be largest. In addition we plotted the 
exact current of the pure system with system size L = I + 1 and a = f3 = 1. 
It seems that the asymptotics for large / are the same as for the bottleneck. 
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Fig. 3. Maximum current in dependence of the bottleneck size for q = 0.6. MC-simu- 
lations are compared with ISA results. The solid line is the current in a homogeneous 
TASEP with L = / + l, a = /3=l and hopping rate q = 0.6. 

In Fig.mthe dependence of the mean density (p) = -^ Yli=i Pi ^^d the current J 
on the entry rate a is plotted for fixed /?. In order to compute these quantities, 
we used the method introduced in [30] that allows to calculate the current and 
densities for an arbitrary set of values of a in one simulation. One observes a 
steep increase of the mean density for the same value of a where the maximum 
current plateau begins. This seems to coincide with the value of a* obtained by 
ISA quite well. The slope increases with system size indicating a discontinuity 
in the thermodynamic limit corresponding to a first order phase transition, in 
contrast to the pure TASEP. The plots also show that at the transition point 
and in the hole PS-phase, the current is maximal. 

As we have argued above, the phase diagram must be symmetric with respect 
to the diagonal a = f3 that yields /?* = a*. With this information we can 
sketch the phase diagram of the TASEP with one bottleneck far from the 
boundaries displayed in Fig. [H Its topology is the same as the one for a single 
defect far from the boundaries [20j , while longer bottlenecks have a larger PS- 
phase than single defects. It also looks similiar to the phase diagram of the 
pure TASEP but here all phase transitions are of first order in contrast to the 
TASEP and the characteristics of the M-phase and PS-phase are different. 

Another procedure to compute the maximum current for a finite bottleneck is 
the finite segment mean field theory (FSMFT) introduced in [2T]. In this ap- 
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Fig. 4. Left: Mean density (p) and current J as function of a. The slope of {p){a) 
becomes steeper for larger systems indicating a discontinuity in the thermodynamic 
limit. The vertical hues show the steepest point of (/>)(«) identified as the phase 
transition a* compared with the ISA results from (fT4l) . 

Right: Schematic phase diagram of the TASEP with one bottleneck far from the 
boundaries. It looks similiar to the phase diagram of the pure TASEP while the 
transition lines to the maximum current (phase separated) phase is shifted to be at 
Of* and /?*. The mean density is discontinuous at these points. 



preach, a segment of ^ sites including the bottleneck of length / < /Ms consid- 
ered. Currents can then be obtained from the eigenvector of the zero-eigenvalue 
of the 2^ X 2Mtransition matrix of this segment [21j. The advantage of this 
method is that the accuracy can be systematically increased by expanding the 
size r of the segment. It can also treat arbitrary combinations of hopping rates 
inside the segment. However, due to the exponentially decreasing size of the 
transition matrix one is currently restricted to segment lengths / < 20. In con- 
trast, the ISA-method, though not asymptotically exact, relies on finding one 
specific root of a polynomial equation of a maximal degree 1 + 2. This makes it 
possible to compute the maximum current for systems with large bottlenecks 
of several hundred sites rather easily. This advantage becomes relevant for dis- 
ordered systems with finite defect site density where bottlenecks of arbitrary 
length can occur. For example the computation of the maximum current of a 
system with 500 consecutive defect sites can be made in less than one second 
on a standarcE] PC to obtain j;^^(/ = 500) = 0.15045 for q = 0.6. Although 
this agrees nicely with the value from MC simulations, J^^^ = 0.15049, it 
is clearly diflFerent from the asymptotic value J^ = q/A = 0.15 for / ^ oc 
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4 Edge effects: One bottleneck near a boundary 



Next we consider a system with a bottleneck near the left boundary i.e. now 
both Li and L2 = / + 1 are of order (9(1) for L ^ oc. The case of a bottleneck 
at distance d from the right boundary i.e. the last slow site is at site L — rf— 1, 
has not to be considered separately since we can deduce the results using the 
particle hole-symmetry 



a ^^ /?, p j ^^ 1 — p j , j ^^ L — j, subsystem 1 ^^ subsystem 3 . (18) 



Since we have only one macroscopic subsystem, namely subsystem 3, the clas- 
sification of phases is slightly different than in Sec. [31 Now the phase of the 
system with bottleneck is basically identically to that of subsystem 3. Phase 
separation can no longer occur since the size of subsystem 1 is microscopic. 
The entry rate a^ of subsystem 3, defined in ISA (see (El)-©), can thus be 
treated as an effective entry rate a^ff := a^ ^ a for the bulk of the system 
which is a homogeneous TASEP. The phase of the full system corresponds to 
the phase of subsystem 3 which we now denote by L', H' or M. 

For Li > 2, we can divide the system into three subsystems in the same 
manner as in Fig. 1. First we consider a system in the low-density phase L'. 
Therefore the current is given by 

J = aeff(l -Q^eff). (19) 



In the steady state the entry and exit rates of a pure TASEP with A^ sites 
are a = j^ and /? = — [13j. The currents in subsystems 1 and 2 have to be 
identical and thus, as in ( ITTl ). 



Jo{a, 1 - p„ Li) = Jo ( ./ , , — , L2] , (20) 

\q{l-Ps) qpf J 

where Jq is the exact current in a pure system given by ( ITOl ) and Ps and pf 
are defined in the same manner as in Sec. [3l Since subsystem 3 is assumed 
to be in the low density phase, the density profile is fiat at its left end, i.e. 
the density is independent of the position if we apply ISA. Therefore we have 
c^eff = P/ = P/+1 ^^d can take J = p/(l — p/) in eq. (I20l) . Furthermore the 
currents in subsystem 1 and subsystem 3 have to be the same which leads to 



Jo{a, 1 - Ps, Li) = pf{l - Pf) . (21) 
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Now we have two ISA-equations, (l20l) and (1211) . as well as the two variables Ps 
and pf = a^ff. For given a, Li and bottleneck length / these equations can be 
solved to obtain the effective entry rate aefr, while the resulting current can 
be obtained from (IT91) . 

For Li = 1 that procedure does not work in the same way because the current 
in a system of size 1 is not defined as one can see in (ITOl) . But even in this 
case, we still have a = J/{1 — pi) which is valid by definition of the boundary 
rates. After inserting J = pf{l — pf)^ this equation together with (I2TI) again 
gives a solvable set of two equations with the two variables Ps and pf. 

For Li = subsystem 1 does not exist and we have only the two subsystems 
2 and 3. Subsystem 2 comprises the sites {1, . . . , L2} and subsystem 3 {L2 + 
1, . . . , L} which includes the bulk of the system. This problem can however be 
solved, by inserting 

a2 = a, /?2 = , as = pf . (22) 

L- Pf 

Analogous to ( JTTl ) we thus obtain the equation 

J 




q' q{l-pf) 



,L2 =p/(l-p/), (23) 



with J = pf{l — Pf). This equation can be solved for the variable pf = aeff 
thus we obtain the effective entry rate and by (IT91) the corresponding current. 

Because of the particle-hole symmetry (IT81) we can transfer these results for 
/?eflp in the high-density phase. The effective exit rate for the bulk /3eff can then 
be determined using (IT9l) - (l23l) . Note that in the low-density phase, /? has no 
infiuence on the bulk of the system, only in a small region near the boundaries. 
The same is valid in a high-density phase for a. 



4.1 Edge effect: discussion 



In Fig. [5] and [6] we have plotted the current in dependence on the position 
of the first defect site for different values of a while (3=1. One observes 
that the position of the bottleneck has a significant infiuence on the current. 
This is called the edge effect^ first observed in [22j. In that work only an 
increase of the current was observed if defects approach the boundary [positive 
edge effect) which can be seen in Fig. [51 In addition, ISA predicts a decrease 
of the current for low entry rates. This negative edge effect is confirmed by 
MC simulations (Fig. [6]). There is also a region of entry rates, where the 
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Fig. 5. Current in a system with a single defect site (left) and a bottleneck of length 
/ = 4 (right) near the left boundary in dependence on its position, for high entry 
rates. The exit rate is /? = 1.0. The system size is L = 500 while the number of 
iterations is 40000000, q = 0.6. The entry rates have been chosen in order to opti- 
mally display the mixed and the positive edge effect. In the system with the longer 
bottleneck (right), the mixed edge effect is almost not visible even for optimized 
parameters. Note that the magnitude is much smaller than in Fig. [6l 
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Fig. 6. Current in a system with a bottleneck near the boundaries in dependence 
of the position of the first defect site, for low entry rates. The exit rate is /? = a, 
system size L = 500 and number of MC steps = 2000000 for bottleneck length / = 1, 
q = 0.6, and a = 0.2 (left), bottleneck length / = 4, g = 0.6 and a = 0.2 (right). 
The scale for the current has been chosen in order to emphasize the form of the 
deviations. 

dependence of the current on the bottleneck position is non-monotonic {mixed 
edge effect^ see Fig. ED ! Simulations for different system sizes indicate that this 
non-monotonic behaviour is not a finite-size effect. Nonetheless, as we can also 
see in the figures, the magnitude of the positive and the mixed edge effect is 
much smaller than the one of the negative edge effect. The ISA results obtained 
from the equations in the last subsection confirm the existence of negative 
and positive edge effect, while the mixed one is not. This indicates that the 
mixed edge effect is caused by correlations at the edge of the bottleneck. In 
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Fig. [6] we plotted the analytical results for comparison. We did not display 
them in Fig. [5] since the deviations due to correlations are larger than the 
positive/mixed edge effect itself. In the appendix it is shown that the negative 
edge effect is predominant in regimes where the current depends significantly 
on the boundary rates. 

Moreover, we see that the current does not attain a plateau value. Instead 
it seems to approach asymptotically the value J = 0.25 (Fig. [Tj). This is 
confirmed by our simulations, where we calculated the current for the very 
high value a = 100 that yielded an effective entry rate of agflf ^0.5. Since we 
assume the effective entry rate to grow monotonic with the real entry rate, we 
can conclude that we have always a^^ < 0.5, thus the maximum current phase 
can never be attained. As we have argued above, there is no phase separation, 
either, if the distance of the bottleneck from the boundary is microscopic, so 
we can state that for a bottleneck near the boundary only the high and low 
density phase can occur. This is plausible since a similar argument for the 
absence of the maximum current phase in subsystem 3 as in Sec. [3] applies. 

Though strictly speaking there is no M/PS-phase if the bottleneck is near a 
boundary, the system still exhibits a kind of crossover. While there is neither a 
plateau region nor a sharp kink in the dependence of the current on the entry 
rate (see Fig. [Tj), for higher entry rates this dependence is rather weak. For a 
bottleneck far from the boundaries we actually have sharp phase transitions, 
so we can say that by approaching the bottleneck to the boundary, the phase 
transition to the M/PS-phase is "softened" into a crossover. 




0.4 0.6 

a 



Fig. 7. Current in dependence on the entry rate a for /3 = 1 and a defect at site 
j = 3. Results are obtained by MC simulations. There is no plateau. 

The transition from low density phase to the high density phase occurs for 
/? = a^ff. This is confirmed in Fig. [8] where we have plotted the average density 
in dependence on a. The jump in the average density marks the transition 
point, which matches quite good the value of a^ff. One observes that the 
transition point is shifted to higher values of a compared to the transition 
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Fig. 8. Average density in dependence on the entry rate for (3 = 0.2 and a defect at 
site 3. The transition to the high density phase is shifted to larger a compared to 
the pure system. 
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Fig. 9. Effective entry rate in dependence on a. Interpreting the y-axis as the beta 
range, this corresponds to the phase diagram with the graph being the transition 
line between H-phase (below) and L-phase (above). 



in the pure system a = (3. In this diagram, a^ff was obtained by the formula 
-. The value calculated by solving the ISA equations (l20l) and (I2TI) 



Q^eff 



-P/+ 



is a little less, but it still yields the correct sign of the shift of the transition 
point. 



Fig. [9] shows the dependence of a^ff on a. Since P = agflf is the transition line 
between high and low density phase, the diagram simultaneously displays the 
phase diagram, interpreting the y-axis as the /?-range. As we see the phase 
transition line calculated in ISA is in good agreement with MC results. 
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5 Two bottlenecks far from the boundaries 



The next step towards disordered systems is the investigation of a system 
with two bottlenecks. This wiU teU us something about the importance of 
"interactions" between the bottlenecks. We distinguish two cases. First we 
will consider situations where both bottlenecks are in the bulk of the system, 
i.e. far away from the boundaries. Then we study edge effects in more detail by 
allowing one bottleneck to be close to one of the boundaries. Unfortunately, 
we cannot use the analytic ISA approach applied in the last sections since 
the equations would be underdetermined. Thus we have to rely on simulation 
results. 



5. 1 Two bottlenecks far from the boundaries 



We simulated systems with two bottlenecks of length /i and I2 with d fast 
sites in between. We focussed on the maximum current phase and determined 
the current J2 • I^ Fig- [13 J2 is plotted as function of the distance d between 
the bottlenecks for different values of /i and I2. One sees that, if the lengths 
of the two bottlenecks differ with h > h^ J2 tends to converge to J*(/i), 
which is the value obtained in a system with only the longer bottleneck. The 
convergence is faster for a larger difference /i — I2 of the bottleneck lengths. 
In this case for a distance of about 5-10 lattice sites, the maximum current is 
almost the same as for a system with only the longer of the two bottlenecks. 
Of course for d ^ the bottlenecks merge, thus we have only one bottleneck 
and ^2(^1^ ^2,0) = J*(/i + h)- If both bottlenecks have equal size, d ^ 00 
converges to the maximum current of a single bottleneck which generalizes 
results of [22] to the case li = I2 > I. 



5.2 Two bottlenecks: Edge effects 



Next we simulated a system with two bottlenecks where one is near the bound- 
ary and one is far away. We concentrated on the case, where the bulk bottle- 
neck is larger than that close to the boundary. 

In Fig. dU we plotted the dependence of the current on the distance of the 
first bottleneck. For comparison we included the results for a single bottleneck 
near the boundary from section [H One observes no significant difference be- 
tween the two datasets. This observation indicates that bottlenecks far from 
the boundary do not have any influence on the current., as long as the current 
is below the maximum current allowed by that bottleneck. 
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Fig. 10. Maximum current in a system with two bottlenecks of length /i (the first 
one) and I2 (the second one) far from the boundaries in dependence on the distance 
of the two defects. The current is determined by MC simulations for a = /? = 0.5. 
In each graph an algebraic fit of the form f{d) = J'^ih) — (d + C)~^ with fit 
parameters C and 7 is included. J*(/i) is obtained from the MC simulations in 
Sec. [3l One observes, that the maximum current converges to the maximum current 
of a system with a single bottleneck, i.e. the longer one. The parameters are /i = 5, 
/2 = 1, g = 0.6 (top left), /i = 5, /2 = 1, g = 0.3 (top right), /^ = 4, /s = 3, g = 0.6 
(bottom left), and /i = 4, /2 = 3, g = 0.3 (bottom right). 



In Fig. [121 however, we see the maximum current in a system with a bottleneck 
of length / far from the boundaries and a defect at site 3. Again, one does 
not see a significant difference: a small bottleneck near the boundary has no 
influence on the transport capacity. 



In agreement with observations already made in [2T] , our results motivate the 
view of a local influence of bottlenecks that yields the possibility to generalize 
concepts of the TASEP with single bottleneck to systems with many bottle- 
necks. Only bottlenecks near the boundaries have influence on the current if 
it is below the transport capacity. We therefore propose that the influence of 
boundary defects can be taken into account through effective boundary rates 
in the same manner as for a single bottleneck near the boundaries. That means 
that below the maximum current we can describe the system as a pure TASEP, 
but with effective boundary rates o^efr,/?eflF depending on the configuration of 
bottlenecks near the boundaries instead of the pure ones. 
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Fig. 11. Dependence of the current on the position of the first defect for a — 0.2 
and (3 — 1.0. Comparison of a system with only one defect and a system with a 
additional bottleneck far from the boundaries. 





' 1 ' 


1 ' 1 ' 1 








+ with defect near boundary 
o withoutdefect near boundary 




S^O.25 

a 0.2 

a 

^0.15 

a 

0.1^ 


1 


^^^^'^^^^^^^^o 


- 


) 5 10 15 20 
length of bottleneck 





Fig. 12. Maximum current in dependence on the bottleneck length for q — 0.6. The 
addition of a defect near the boundaries does not alter the maximum current. 

6 Discussion and Outlook 



We have investigated the TASEP in the presence of a bottleneck of / con- 
secutive slow sites. Apart from computer simulations we have developed an 
approximate analytical approach called interacting subsystem approximation 
(ISA) that makes use of the exact results obtained for the pure TASEP. ISA 
turns out to be an efficient and accurate method for systems with a single 
bottleneck. It can be applied to other systems as long as exact results are 
available. 

Therefore we have divided the system into subsystems consisting of fast and 
slow sites only, where the exact solution p!3l[T4] for the homogeneous system 
applies. In the treatment of the coupling of these systems certain correla- 
tions are neglected. Nevertheless the predictions of the analytical approach 
are in good agreement with the simulations, while they slightly underesti- 
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mate the real current. Our method yields much better results than e.g. the 
mean field approach in [20j for single defects. For longer bottlenecks it is much 
more computationally efficient than the approach by Chou and Lakatos [21], 
which produces results of a similar accuracy. The finite segment mean field 
theory (FSMFT) introduced in [21 J takes into account only particle correla- 
tions within the bottleneck section and a certain number of sites adjacent to 
it. This approach is more ffexible in the kind of systems that can be treated, 
but it is algorithmically more complex and less efficient since it requires the 
determination of the leading eigenvector of a matrix that grows exponentially 
in the segment length /. As a consequence it can not be applied to systems 
with longer bottlenecks (/ > 20). In contrast, the ISA method relies on finding 
one root of a polynomial whose degree is approximately /. This root can be 
isolated by physical constraints. For this problem fast algorithms exist mak- 
ing it possible to compute the current for very large bottlenecks. Moreover the 
accuracy of the ISA even increases for longer bottlenecks. Since later [2lj we 
want to use our results for disordered systems with finite density of defects, a 
method that can also be applied to long bottlenecks is required and the ISA 
appears to be a suitable procedure for this task. 

The ISA can be applied for calculating the transport capacity (maximum 
current) and the transition point to the phase separated phase of a system with 
a bottleneck far from the boundaries, as well as the current of a system with 
a bottleneck near the boundaries. If the bottleneck is far from the boundaries, 
the main effect on the current is a lowering of the maximum current. It is 
independent of the bottleneck position, but things change if the bottleneck is 
located near one of the boundaries. Here interactions with the particle input 
or output become important [edge effect^ for single defect sites see [22j). This 
leads to a change of the current. ISA reproduces the increase of the transport 
capacity if defects are moved towards the boundary which was found by Dong 
et al. [22]). Additional to this positive edge effect, ISA predicts that for lower 
entry rates the current is decreased by a bottleneck near a boundary compared 
to the situation where it is deep in the bulk. This negative edge effect is 
confirmed by numerical simulations. The simulations also show that there is a 
regime where even a non-monotonic dependence on the bottleneck position is 
possible. This effect, however, is not reproduced by ISA which indicates that 
it has its origin in correlations that are neglected in the approximation. 

A bottleneck near the boundary also has crucial effects on the phase diagram. 
Indeed it leads to an absence of a maximum current-/phase separated phase! 
One can say, the bottleneck near the boundary softens the phase transition 
into a crossover. The remaining phase transition line between low and high 
density phase, on the other hand, is distorted compared to the pure system. 
This observation motivates the concept of effective boundary rates: the effects 
of the boundary bottleneck are encompassed as effective boundary rates and 
the system can be described as a pure one with the effective rates instead of 
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the real ones. One can calculate the relation between effective and real rates 
by ISA which then can be used to determine the phase transition line. 

In future work [2l] we will use the present results to derive an effective theory 
for disordered systems which have a finite density of defect sites even in the 
thermodynamic limit. The results obtained here for systems with two bot- 
tlenecks indicate that an accurate description of the disordered case can be 
obtained, since the "interaction" of bottlenecks is rather local which is in agree- 
ment with observations in [2T] and [22]. For large distances of bottlenecks, the 



maximum current in the TASEP is dominated by the longest bottleneck while 
(smaller) bottlenecks even close to the boundaries do not have any influence 
on it. This observation supports the conjecture made in [261127] for periodic 
systems with quenched disorder that the longest bottleneck is dominating the 
transport capacity, and hence the transition to phase separation. Therefore 
it will be important to have an analytical approach that can deal efficiently 
with very long bottlenecks. In contrast, if the system is not in the maximum 
current phase, we have found that even large bottlenecks in the bulk do not 
have an influence on the current. These observations indicate that we can treat 
boundary effects and transport capacity independently for systems with many 
bottlenecks by using the concept of effective boundary rates and the ISA also 
in this case. 

Natural applications of the results presented here are to biological transport. 
The TASEP is the basis for the description of protein synthesis and also ac- 
tive intracellular transport of motor proteins along microtubule In both cases 
disorder plays an importat role through the inhomogeneity of the mRNA tem- 
plate [331135] or microtubuli associated proteins (MAPs) that attach to micro- 
tubuli and can form obstacles for the motion of mobile motor proteins like 
kinesin 



Appendix: More on edge effects 



We now want to determine a condition that the current decreases for a bottle- 
neck approaching a boundary {negative edge effect). Again, due to the particle- 
hole symmetry (IT81) it is sufficient to consider a bottleneck near the left bound- 
ary. We have to distinguish if the current for a defect at a small distance 
Li = 0{1) from the boundary is larger or less than the current for a defect far 
from the boundary, Li ^ oc. Though the dependence of the current on the 
position of the first defect is not always monotonic, for an approximation it 
is sufficient to compare the current J{Li) for one position of the defect near 
the boundary with the current J(oc) for a defect far from the boundary, since 
the mixed edge effect (see Sec. Hj) is quite weak and only in a small parameter 
region. For simplicity we take Li = 2. 
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Using (HD and © yields 

Pi + api + a^ + api + a^p\ 



where /?i is the virtual exit rate of subsystem 1. The edge effect is negative if 

J(a,A,2) < J(a,/3i,(X)). (25) 

In the L-phasO we have J{a^ Pi^ oc) = a{l — a) and f l25f l is fuUfiUed if 

/?i < 1 - a . (26) 



Though we cannot determine /?i exactly, we can state that /?i is less than in 
a pure system with q = 1. However, in a pure system there is a flat density 
proflle without correlations, and we have exactly /?i = 1 — a, so f l26l ) is always 
fulfilled. Therefore in the L-phase and, due to particle-hole-symmetry, in the 
H-phase, the edge effect is always negative. 

In systems with many randomly distributed bottlenecks it is very unlikely that 
longest one is near the boundary. Thus the transition to the PS-phase, where 
the current is independent of the boundary rates, will occur for lower rates of 
a(/3), where the negative edge effect is predominant. 
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